Breaking of general rotational symmetries by multi-dimensional classical ratchets 
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We demonstrate that a particle driven by a set of spatially uncorrelated, independent colored 
noise forces in a bounded, multidimensional potential exhibits rotations that are independent of 
the initial conditions. We calculate the particle currents in terms of the noise statistics and the 
potential asymmetries by deriving an n-dimensional Fokker-Planck equation in the small correlation 
time limit. We analyze a variety of flow patterns for various potential structures, generating various 
combinations of laminar and rotational flows. 
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I. INTRODUCTION 

One of the corner-stones of equilibrium statistical me- 
chanics is the second law of thermodynamics, which pre- 
cludes extraction of pure work out of a heat source (e.g. a 
thermal noise) without an accompanying change of state 
0. Fluctuation-dissipation theorem makes it quite im- 
possible for a system to extract work out of a noise source 
in equilibrium with it, even if the system is in a potential 
with a built-in directionality. The impossibility of such a 
thermal noise-induced rectification was lucidly explained 
by Feynman |2[] using the concept of a 'ratchet', a device 
whose static potential is periodic, but with an asymme- 
try within each period (such as a sawtooth potential). 
Although motion in a ratchet is easier in one direction 
than the other, it is impossible to exploit this asymmetry 
to drive a particle in the potential using a simple ther- 
mal environment. This is because the probability of a 
noise-induced jump over a barrier in a ratchet depends 
only on the barrier height, and is therefore the same to 
the right and the left, irrespective of the different slopes 
in the two directions. 

Directed motion in a ratchet potential requires an ex- 
ternal source of energy that is out of equilibrium with 
the system, thus negating the necessity to obey the 
fluctuation-dissipation theorem 0,0. A trivial exam- 
ple is a fully correlated noise source such as a determin- 
istic unidirectional force exerted to an axle attached to 
a ratchet-wheel. Remarkably, the rectification persists 
even for a non-equilibrium noise with zero time-average. 
The minimal condition for a ratchet to operate is broken 
detailed-balance (such as a "colored" or temporally corre- 
lated external noise) which in conjunction with the spatial 
asymmetry within each ratchet period, and broken time- 
reversal symmetry (dissipation) leads to directed motion. 

The successful extraction of work out of a non- 
equilibrium source of energy has far-reaching implica- 
tions. Thermal ratchets are not limited by energetic re- 
strictions associated with equilibrium statistical mechan- 
ical principles 0, IE HI • Massive (underdamped) ratch- 
ets exhibit a parametric current reversal that could be 
useful for continuum mass separation [Toj and designing 



'molecular shuttles' [Tl|. Furthermore, ratchet motion 
is considered to be a possible explanation for the long- 
range cellular transport of motor proteins [T^|. On the 
experimental front, Brownian ratchets have been demon- 
strated in the rectified motion of polystyrene spheres and 
a drop of mercury [l3| , and in current rectification in a 
DC SQUID 0. 

A ratchet-induced rectification using a two or three di- 
mensional apparatus has been investigated recently, pri- 
marily with the aim of separati ng p articles of different 
diffusion constants or sizes ITol Il5l. Though the par- 
ticle motion in such studies [lfi Il5l Hif is in multi- 
dimensions the ratchet effect is usually along a single 
direction. Relatively little investigation has been done of 
systems where the ratchet effect exists in multiple spatial 
directions \vA fl8| . Extending ratchet motion to higher 
dimensions involves more than just a simple extension 
of one-dimensional arguments. One can envisage break- 
ing of higher symmetries such as rotational symmetry, by 
the broken detailed balance, leading to a rich structure of 
loops and vortices. Generalizing the concept of a 'coordi- 
nate', one can map the motion of a stochastic particle in 
phase-space (semiclassical evolution in a multi-quantum 
well system for example), or in chemical coordinates, into 
a n-dimensional motion in real space. The interplay of 
rectification and dimensionality can lead to very interest- 
ing flow patterns. For instance, a judicious combination 
of one-dimensional ratchets can lead to steady-state rota- 
tions, even if the noise sources themselves are restricted 
only to apply along the orthogonal directions and are un- 
correlated with each other. A simple realization in two 
dimensions (2D) is as follows: Consider Fig. ^where a 
two-dimensional rotation due to ratchet motion is shown 
schematically. In the presence of spatial asymmetry in 
one-dimension, a time-correlated noise is known to pro- 
duce a drift 0, EI; the direction of the drift is deter- 
mined by the sense of the potential asymmetry. For a 
potential in multiple dimensions, the sense of the po- 
tential asymmetry along one coordinate can be reversed 
by varying the other coordinates, leading to a change in 
sign (or reversal of asymmetry) of the potential. As in- 
dicated in Fig. ^ the coordinate-dependent reversal of 
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x-drift 




FIG. 1: Schematic description of rotation over one unit cell of 
a two-dimensional periodic potential, caused by spatial asym- 
metry and temporal correlations. One-dimensional drifts are 
produced by asymmetric potentials in x and y, in conjunction 
with time-correlated noise. The drifts along x and y direc- 
tions switch signs owing to coordinate-dependent changes in 
overall sign of the potential asymmetry, and in combination 
produce rotation. A specific example of a potential V(x,y) 
that generates these flows is shown in color. 



the one-dimensional drifts could then conspire together 
to generate a steady-state rotation. Contrary to rota- 
tion generated by purely potential forces, determined by 
the initial conditions of the particle, the sense of our ro- 
tation is independent of initial conditions and is given 
by the combination of the potential asymmetry and the 
noise statistics. Furthermore, removing either the poten- 
tial asymmetry or the correlation in the noise destroys 
the rotation. In effect, we have thus produced a rota- 
tional motion along a cyclical track using combinations 
of one-dimensional ratchets, although the x and y noise 
forces are totally uncorrelated. A generalization of this 
process can generate rotations for any asymmetric po- 
tential in x and y that is non-separable in the individual 
coordinates. Rotation then becomes a necessary outcome 
of motion in such a potential • 

In this paper, we develop a formalism for dealing with 
classical ratchet motion in n-dimensions. We generate 
a bona-fide Fokker-Planck equation (FPE) in terms of a 
systematic expansion in the correlation time for a Gaus- 
sian distributed noise. The work is presented as follows. 
Section II sets up the functional calculus in arbitrary di- 
mensions that allows us to derive a bona-fide FPE in the 
presence of colored noise. The formalism in itself is non- 
perturbative in the noise color. In section III we derive 
the n-dimensional FPE. In section IV we specialize to 
the case of one-dimension and analyze the ratchet mo- 
tion. We also discuss the various length and time scales 
built into the dynamical process of a particle in a poten- 
tial driven by colored noise. In section V we present a 
time-dependent solution to the n-dimensional FPE, just 
before steady state is reached. In section VI we analyze 
the steady state part of the solution of section V and de- 



scribe the variety of flow patterns that may arise in two 
and three dimensions. We conclude in section VII. Ap- 
pendix A gives some algebraic details of the calculations. 
Appendix B discusses questions about convergence. 



II. FUNCTIONAL CALCULUS FOR THE 
PROBABILITY IN MULTIPLE-DIMENSIONS 

Consider an overdamped classical particle in a multi- 
dimensional potential U(r) driven by an external time- 
correlated noise with force jf{t), where 7 is the damp- 
ing constant. The motion of the particle is governed by 
the Langevin equation, which is essentially Newton's law 
without the inertial term: 



r(t) = W(f) + f(t) 



(1) 



where r is the position vector of the particle, jW = 
—VU if) is the force exerted by the potential on the par- 
ticle, and the dot denotes a derivative with respect to 
time. The stochastic process described above depends 
crucially on the statistics of the noise /. We assume that 
the noise has a Gaussian probability distribution with a 
correlation time t£ and strength D t along the zth coor- 
dinate direction. Furthermore, we assume that the noise 
along any two orthogonal coordinate directions are un- 
correlated. The probability distribution P(r,t) j2^| can 
be described in terms of a functional integral over differ- 
ent realizations of / as follows: 



P(f,t) = / VfP[f\8{r-r {€))., 



P[f] = Vexp 



\J J dsds'Y^Kij(s- s')fi(s)fj(s') 
</i) = 0; W)fjV)) = fe^W 

\ 'i J 



(2) 



The Dirac delta function in Eq. equates the arbitrary 
position variable r with the functional form r{t) stipu- 
lated by Eq. . The dimensionless correlation function 
Cj can in principle have multiple time-scales built into 
it. The formalism we develop can be easily generalized 
to take such extensions into account. 

We generalize the functional calculus outlined by Fox 
[T^ | to develop a bona-fide Fokker-Planck equation for the 
probability distribution P(r,t). Using Eqs. QJ and @ 
we get 



dP(r, t) 
dt~ 



-JvfPlfi^l^ 



UVi 



(3) 



where Q t (r) = / VfP{f}5(r-r(t))f t (t). Using results de- 
rived in Appendix A, we can write the product P[f]fi(t) 
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as follows: 



P[f\h{t) 



D, 



i 



dsPlfjf^SuSit - s) 



i 



-i / dsP[/]/,(s) / ds' K u (s - s')a, 



Substituting in Eq. |3] 
dP(r, t) 



where 



dt 



i(t) 



= -v • (WP) - 



d 2 Q tJ 
dridrj 



(4) 



(5) 



x VfP[f]S(r-r(t)) 



5rj(t) 



(6) 



At this stage, we see that in the presence of color, the 
stochastic process is non-Markovian, which means that 
the value of P at time t depends on earlier instants of time 
s' through the kernel K{t — s'), or equivalently, through 
the correlation function C which represents the inverse 
of K (Appendix A). To evaluate the role of color, we now 
need to evaluate the functional derivative in Eq. [B] This 
is best done with the help of the Langevin Eq. ^ for fj : 



S rj (t) 



6+i (t) 
Sfi(s') 



dt [dfi(s')_ 

= mo lWl+m 



k 



where we define Mi 
above equation is: 



Sfi(s>) 



dWi(t)/d rj 



(7) 

The solution to the 



Shis') 



H(t - s')SijT exp 



dsM(r(s)) 



(8) 



j ' 



as can be verified by direct substitution. In the above, T 
is the time-ordering operator and H(x) is the Heaviside 
step-function (defined to be unity for positive x and zero 
otherwise) . Substituting this functional derivative in Eq. 

we manage to reduce the equation for P(f,t) into 
the form of a continuity equation: 

P = — V • J, with 
d 



Ji 



mp- w [e iP] , 



d r°° t' ( \ 

6, = ^/ dt'Ci(—)T exp / dsM{ 
T t Jo T t \ Jt- 1' 



(f(s)) 



(9) 



III. GETTING A BONA-FIDE FOKKER-PLANCK 
EQUATION 

To get the steady state (t — > oo) limit of the above 
equation, we expand the integral inside the exponential: 



rt f>2 f) 2 W ■ 



fk] + 0{t' 2 ) 



(10) 



The expansion does not introduce any singularities, as 
discussed in the Appendix B. After changing variables 
t' — > x = t' /t£, we can now rewrite 0i as: 

/>oo 

®i=Di / dxCi(x)T(exp [xt^M -x 2 {t^) 2 R/2 
Jo 

+ 0(x 3 (T?f)}) ii (ll) 

where the matrix R has components i?y given by Rij = 
J2k Wkd 2 Wi/drjdrk- The components of M and R have 
dimensions of w Uq/jL 2 (Uq is the maximum height 
of the potential and Li is the length scale of variation of 
the potential in the i-th direction, which equals its pe- 
riod for a periodic potential) while x is dimcnsionless. For 
small correlation times rf <C (time scales described in 
detail in the next section), we can further Taylor expand 
the exponent to give us the effective diffusion constant 
Qi in terms of the noise statistics defined by {m^}, where 
fi l n = L dxCi(x)x n is the nth moment of the noise corre- 
lation function. For well-behaved functions with rapidly 
vanishing higher moments, terms like ir n (7f )" can be ig- 
nored for small r? and large n (cf. Appendix B). Then 
we can truncate the equation of motion for P(r, t) to sec- 
ond order, leading thereby to a bona-fide Fokker-Planck 
equation. For small correlation time r 2 c , the equation 
reads: 



dP(f, t) 



dt 



E 



e, 



M, 



Ri 



D; 



dri ^—^ 9r\ 

(rf) 2 



1 + (j\T?Mi, 



lA{R-M 2 ) 



= dWjjt) 



E^ 



d 2 W t 
drjdru 



(12) 



The above set of equations are the central equations 
for all our analyses. The effect of the noise correlation 
shows up in the effective diffusion constant Oi ; which 
picks up a position dependence J[2$J in a well-defined man- 
ner through the potential gradient terms M and R. The 
statistics of the noise shows up through the moments 
{/J- l n } of the temporal correlation function. 
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IV. APPLICATION: 1-D RATCHET 

Having established an approximate though bona-fide 
Fokker-Planck equation for arbitrary correlation func- 
tions d and arbitrary dimensions at small correlation 
times, we now use our results to obtain the dynamics 
of a classical particle in various kinds of ratchet poten- 
tials. As a first step, we calculate the steady-state cur- 
rent density in a one-dimensional periodic potential using 
our Fokker-Planck formalism for exponential correlation 
(fii = 1, /i2 = 2). In one-dimension, 8 can be written as: 



n 



L dx^M- [1 - t c W'(x) + (t c ) 2 W(x)W"(x) 



D 



+0((r c ) 3 )] 



(18) 



The first and second terms in the integrand are propor- 
tional to exact derivatives (of U and W 2 /2 respectively), 
and vanish due to periodic boundary conditions, so that 
J oc (t c ) 2 . This allows us to replace 1 — exp (—cj)(L)) in 
Eq. (|16fl by the small quantity <fi(L). The current density 
is then given, to leading order in t c , by: 



6 = D [l + t c W - (t c ) 2 {WW" - W' 2 }] 



(13) 



to second order in r c and where prime denotes a deriva- 
tive. The expression for to first order in r c is well 
known in the literature [l9| . However, we need to retain 
the second order corrections in r c to get any nontrivial 
current density out of the noise, as we shall shortly see. 

At steady state (P = 0), the Fokker-Planck equation 
reads: 



d 4- = o 

ax 



J 



WP - — [6P] 
ax 



(14) 



The continuity equation dictates a constant (x- 
independent) current. Using an integrating factor 
exp {—<t>) where <j)(x) — J Q dyW(y)/Q(y), '0' being an 
arbitrary point on the x axis, we solve the above first 
order ordinary differential equation for P: 



P{x)Q{x)e-' l3{x) = P(0)6(0) — J I dye 



-<P(y) 



(15) 



Imposing periodic boundary conditions -P(O) = P(L), 
0(0) = <d(L), one then gets the following equation for 
the drift current J: 



J = P(0)9(0) 



1 - e 



I I dxe-^ x) (16) 
o 



We can solve for P(O)0(O) by normalizing P{x) within a 
period with a normalization constant N representing the 
number of particles within a period at steady-state. This 
yields finally: 



J 



N [l 



-<j>{L) 



4>{L) 



[W(x)/Q{x))dx 



J = (r c f 



N Jq dxW 2 (x)W"(x) 
dxe- u ^)hD ft faeUWhD 



(19) 



A non-zero drift current density in a ratchet is gener- 
ated because the following two quantities don't vanish: 
(a) correlation time (r c ^ 0), which signifies a non- 
equilibrium noise that breaks detailed balance, and (b) 
potential asymmetry , meaning that there exists no Ax 
such that U(x + Ax) — —U(x). This makes the inte- 
gral in the numerator of Eq. 1191 non-zero although the 
integrand itself is periodic. 

The dynamics of the particle has three time-scales built 
into it: 

(i) Correlation time, t c . This is the time governing the 
rate of loss of 'memory' in the noise. A convenient way 
to view the correlation time (as we establish later in this 
section) is an effective partitioning of the dynamics such 
that for times less than the correlation time the motion 
is ballistic, while for times greater than the correlation 
time the motion is diffusive. 

(ii) Diffusion time, t d — L 2 /D. L is the typical length 
scale built into the potential. For a periodic potential, 
for example, L denotes the period. t d describes the time 
taken to diffuse over one unit length scale L of the po- 
tential. 

(Hi) Drift time, r 7 = 7L 2 /Uq. Uq is the maximum height 
of the potential, which serves as a typical energy scale in 
the problem. The drift time describes the amount of time 
for an overdamped particle to drift from the maximum 
to the minimum of the potential over a distance ~ L, 
and serves as the time-scale for noise-free motion in the 
system. 

Let us calculate J for a 1-D periodic asymmetric po- 
tential with period L and height Uq, and impose periodic 
boundary conditions P(0) = P{L). Using Eq. ^J] for J 
we get the 1-D drift current as: 



dx €§) K -I 1 - ^ (L) ] Jo* d y^ (v) } 

. . (17) 

For small correlation times r c , using the expression for 
<d(x) from Eq. I|13|) . one can expand <j>(L) as follows: 



j (r c r D ) 2 

(tT) 5 9 \ r " 



(20) 



where g{x) is a factor that depends on the geometry of 
the potential V(x). Fig.|3shows the steady-state current 
density for a specific example of a periodic potential for 
varying correlation times r c . The current density is posi- 
tive, tending to drive the particle out of the potential well 
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FIG. 2: Potential U(x) = Ax 13 (l-xY (solid line) periodically 
extended, and the corresponding steady-state probability dis- 
tributions for t c = 1 x 10~ 5 (dash-dot), 1 xlO -3 (dotted) 
and 1 xlO -2 (dashed). Here the potential height Uo = 1 and 
D — 1. As t c increases, the probability peaks to the left in 
the right well, indicating a positive current density, in the 
direction of lower restoring force. 

in the direction where the restoring force —dU{x)/dx is 
less. 

A heuristic argument can illustrate the origin of the 
directionality of the current. For a correlation function 
(/(t)/(0)) = (D/T c )C(\t\/T c ), one can use Langevin's 
equation to get: 

(r 2 (t)) = 2Dt / dYC(Y) -2Dt c / dYYC{Y) 
Jo Jo 

= 2DC(0)t 2 /r (t <r c ), 

= 2D[t-mT] (i>r c ) (21) 

where fi± is the first moment of the noise. The above 
scaling of the particle coordinate with time means that 
the particle motion can be looked upon as ballistic for 
t <C t c and diffusive for t ^> t c . In absence of correlations 
(t c = 0), the particle jumps over the barriers purely due 
to diffusive transport, and this has equal likelihood in 
either direction, since the barrier height is the same to 
the right and the left. However, for a correlated noise, 
there is a net drift over a time r c that pushes the particle 
ballistically more to the right than the left, since the 
restoring force is less to the right (see Fig|2Jl. At the end 
of the drift process therefore, the particle has reached a 
higher elevation to the right than the left. Thereafter, 
the barrier height is smaller to the right, so there is a 
higher probability of crossing it to the right. Thus the 
net drift in the positive x direction is a consequence of 
drift assisted thermal activation, produced by the presence 
of correlation in the noise. 
Fig. |21 shows the steady-state probability distribution 




Barrier Height (arb. units) 



FIG. 3: Steady-state current density for a U(x) = 
Ax 13 (1 — x) 7 , plotted for varying values of the maximum po- 
tential energy Uo. The peak structure arises out of an inter- 
play between increasing ratchet effect and increasing backflow 
of current as the potential maximum is increased. As Uo in- 
creases from 0, initially the ratchet effect dominates and when 
it becomes too high the current back-flow dominates giving 
rise to a peak in the ratchet current at intermediate values. 



function for varying values of the correlation time r c . For 
white noise (r c = 0), the probability distribution is given 
by the Maxwell-Boltzmann result P oc exp [—U(x)/ r yD]. 
As the correlation time is cranked up, however, the prob- 
ability distribution within each well of the potential pro- 
gressively shifts in weight to the left. This corresponds 
to a non-vanishing steady-state current density J in the 
positive direction. The direction of the current density 
can be obtained easily from Eq. {H}, whereby the peak 
(dP/dx = 0) of the probability density is given for small 
correlation times (0 rj D) by the point where J = WP. 
Since P > 0, therefore for positive current J > 0, we have 
—W = dll/dx < 0, meaning the peak shifts to the left. 
Note that the directionality of the ratchet motion de- 
pends on the specific nature of the correlation, and can 
be different depending on whether the time-correlation 
is incorporated into the driving force, or into a fluctuat- 
ing version of the periodic asymmetric potential we have 
been dealing with. 

The current density plotted vs. the barrier height 
shows a maximum, as is seen in Fig. 3. For small bar- 
rier heights, there is an equal chance of jumping over the 
barrier to the right and the left, so the current is small. 
Similarly at large barrier heights, the chance is equally 
small, so the current is small once again. Thus the peak 
in the current arises out of a competition between ratchet 
motion to the right, and current back-flow to the left. 

One dimensional ratchet motion has been observed for 
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polystyrene beads in an aqueous solution driven by a si- 
nusoidal current between lithographically patterned elec- 
trodes [Hjj . In addition, ratchet motion has also been 
observed for beads in a "flashing potential" generated by 
a laser system modulated suitably by a chopper [3J. In 
the flashing geometry, the beads are trapped in deep po- 
tentials while the laser is on, and allowed to diffuse sym- 
metrically while the laser is turned off. On subsequently 
turning on the laser, the beads slide back into the wells 
of the pieces of the potential whose basin of attraction 
they are individually in. For an asymmetric potential, 
the beads distribute asymmetrically. The entire process 
of turning on, switching off and turning on again leads 
in effect to a net drift of the particles out of the well in 
the direction of the steeper slope (opposite to the drift 
direction for the noise force-driven ratchet). Such drifts 
were reported by direct visual observation. 

Our set-ups can be imagined to be the same, except 
that the driving force is neither an AC signal nor a 
switching potential. Instead, we imagine a static asym- 
metric potential generated as above, but driven by a col- 
ored noise generator (the simplest example is a white 
noise passed through an RC low-pass filter). Imagine a 
potential with period ~ I^m (a = 0.7/im, b — 0.3/im) 
and fluorescent charged polystyrene beads (~ 0.07 - I 
/im) diameter in an aqueous solution at room tempera- 
ture. The energy of the potential Uq is set to ~ 75 meV. 
For this set of parameters, t d « 3 s and t 7 a 1 s. For a 
piecewise linear potential, /(t d /t 7 ) « t d /t 7 . Assume a 
bandwidth of the colored noise of ~ 40 Hz (these param- 
eters satisfy t c <§; r 7 < t d ). Then t « 10 hours. Such 
a slow drift of the fluorescent beads under colored noise 
should be readily observable with a microscope. Subse- 
quently, measuring the current density for varying values 
of Uq should generate a graph similar to Fig. |3 the posi- 
tion of the peak depending on the shape of the potential, 
but the tail varying as exp(— Uq/jD). 

Let us calculate the efficiency of the ratchet process. 
Essentially, the energy for the directed motion is ex- 
tracted from the correlated pieces of the noise force. The 
input power, given by the product of the force of the noise 
and the velocity obtained in the absence of the potential, 
is given by: 



p 

- 1 1 t 



(7/ x /} 



(22) 



To calculate the output power, we use the generated drift 
velocity, and the force —dll/dx(x) required to overcome 
a potential barrier and get: 



Pout 

This yields an efficiency 



T U » 



(23) 



V 



Pout (t ) ( T ) , D I 1\ /nA\ 



For the typical values cited above for our proposed exper- 
iment, this corresponds to rj « 10~ 5 . The transduction 
is not very efficient, although there have been propos- 
als for ratchet mechanisms that perform almost with 100 
percent Carnot efficiency poj. 



V. TIME-DEPENDENT PROBLEM 

The steady-state probability distribution obtained in 
Eq. El is in effect the long-term (t — > oo) solution to 
the full time-dependent Langevin problem. To solve for 
the transient response P(r, i), let us expand the spa- 
tially periodic variables W and O; in their Fourier modes: 



ik ■ 



Expanding 



P(f, t) = J2k Pk elk ' f %W%t i we g et from E q- 



J2 - iuj k p k e 



ik ■ r — i^jX 



k.k' 



k 



(25) 



Changing summation variables on the right from k, k' to 
k, k + k! , we have: 



itort 
e k x 



ik ■ r 



ik' ■ 



Setting the terms in square brackets equal to zero, and 
performing an inverse Fourier transform, we get the dis- 
persion relation: 



-fc 2 (0) - ik ■ (W), 



(27) 



(r 7 ) 



where (. . .} denotes a spatial average over the period L. 
Substituting back into the definition of P, and including 
the steady-state solution as the t — > oo, i.e., the long 
wavelength (k = 0) limit, we get: 

P{r,t) = ^P^expjifc • [f- (W)t] - k 2 (Q)t} +P st (r), 

fe#0 

(28) 

as one expects for a drift-diffusion equation. Since the 
average effective diffusion constant (O) is positive, the 
probability distribution decays with time to the stable 
steady-state solution. 

The "initial conditions" for the Fourier coefficients Pg 
in Eq. |2H|are set by the value of P(r, to), where to 3> t c . 

WegetP£ = / d 3 rP(r,t )exp -k-r + k 2 (B)to /L 3 . At 

times less than r c , memory effects become important, 



7 



the temporal evolution is non-Markovian and the corre- 
sponding equation for the probability distribution P{x, t) 
cannot be truncated to second derivatives in x to con- 
struct the bona-fide FPE. This means that for small 
times the transient equation is not of a drift-diffusion 
form, but depends in fact on higher correlations. This 
disallows the use of the FPE structure to start from an 
initial condition at t — and propagate to steady-state. 
However using an intermediate distribution at to (with 
t 3> t c ) is allowed. In order to obtain this intermedi- 
ate distribution from an initial condition at t = 0, one 
actually needs to work with Eq. [5J which has a more 
complicated time-dependence and associated memory ef- 
fects built into it. However, for times much larger than 
the lifetime of the memory effects (~ r c ), Ea.1281 should 
work well. For times larger than the correlation time t c , 
the approach to equilibrium is governed by the diffusion 
constant averaged over a potential period, (0) , while 
the associated drift is governed by the corresponding av- 
erage potential force, (W), which is zero for a periodic 
potential. 



VI. MULTIDIMENSIONAL ANALYSIS; 
ROTATIONS AND PATTERNS 



Let us generalize our results of the previous sections 
to higher dimensions. We consider a periodic potential 
so that U{r)\ ri —o = U{r)\ ri= i Ji , where L, is the period 
of the potential along the ith direction. Imposing pe- 
riodic boundary conditions along the ith direction and 
integrating Eq. ^] leads to an integral equation for J,: 



dnJi^e-M^ = [P(r)9 4 (r)] 



1-e" 



Ti=Li 



where 



<kif) 



dziiWiW/eiiz)}, 



(29) 



(30) 



In one-dimension, steady-state implies constant current 
density which allows us to pull J out of the integrals 
(as shown in Eg. I15f) and solve for it, with the boundary 
value of P{r) at = and Li being fixed by normaliza- 
tion. The situation is quite different in multi-dimensions. 
At steady-state, the current density is not a constant in 
general. However, we can still make a few observations 
that lead to nontrivial conclusions: (i) the right hand 
side of Eq. I2H1 is m general not identically zero, (ii) the 
integrand on the left hand side is a product of Ji{r) and 
a positive definite quantity, and (iii) definition of steady 
state (P(r, t) = 0) implies V • J(r) = 0. The first two ob- 
servations imply that J cannot be identically zero every- 
where. Combining this with the third observation leads 



to the unavoidable conclusion that V x J is not identi- 
cally zero over one period of the potential (excluding the 
trivial case where J is a constant vector). In other words, 
there necessarily are local rotational patterns. 

The necessity of color and potential asymmetry in our 
arguments is now easily seen. Analogous to Ea.lTSlof the 
one-dimensional case, for small correlation times rf, we 
get: 



f L i W.(Z) 

<>AL,)= / dzi-j^- [l-^M u (z)+ 



+0((7f) 3 )]. (31) 

Terms of zeroth and first order in t£ in Eq. 1311 are zero 
since the integrands in them are proportional to exact 
derivatives (of U and Wf respectively), and vanish due 
to periodic boundary conditions. Hence, 4>i{Li) oc (t?) 2 . 
From Eq.|^ we see that 4>i{Li) = either for white noise 
(i.e. when rf = 0) or when the potential is symmetric 
within a single period (i.e. the net integral multiplying 
( T i) 2 vanishes). If either of this holds then our obser- 
vation (i) is invalidated leaving open the possibility that 
Ji{r) = everywhere. 

Having established the existence of rotations, we need 
to solve Eq. 1121 numerically in n-dimensions with given 
boundary conditions to get specific flow patterns for J(r). 
For the purpose of illustrations, we specifically adopt 
the following simplifications: (a) we concentrate on two- 
dimensions (i — x,y), where Eq. 1311 becomes: 



4>x{Lx,y) 



(r c x f 



D , 



3_m| 
4 



2 



dxW x (x,y) 
d 2 W x 



dx 2 



2 

4-w, 



dW x 
dy 
d 2 W. 



v dxdy 



;(32) 



(b) Next we will consider the specific case of exponen- 
tial correlation function (n\ = 1, n\ = 2); (c) finally, 
we consider a restricted class of potentials and noise 
such that U[x,y) = U(y,x), D x = D y and t^. = rf. 
At steady-state, this implies that J x (x,y) — J x (y) and 
J y {x,y) = J y (x), i.e. the x component of vector J(x,y) 
is only a function of y and its y component is only a 
function of x. Before proceeding we emphasize that our 
arguments (following Eq. I29|) showing the necessary exis- 
tence of rotations do not depend on this restricted class 
of noises and potentials; the class of potentials is adopted 
just to simplify the algebra for illustrative purposes. The 
calculation of J x {y) and J y (x) can now proceed smoothly 
by observing that J x (or J y ) may be pulled out of the in- 
tegral in Eq. [2H The total current density J{r) is then 
obtained by solving the set of Eqs. - The 

current density J x (y) depends on the Dirichlet bound- 
ary conditions P(0,y), which we will set to constant 
(P(0,y) = P(x,0) — constant), since it allows us to get 
simple flow patterns. 
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Flow pattern Coupling mechanism 

Rotation Ratchets coupled in x and y 

Laminar flow Decoupled ratchets in x and y 

Rotation + net drift Coupled ratchets asymmetric under 

x <-» -x, y <-> -y 




0.2 0.4 0.6 0.8 1 

x 

FIG. 4: Contour plots of the potential U(a;, y) — sin(x)sin(y) 
—asin(2x)sin(2y) with a = 1, and all lengths expressed in 
units of 2-7T. White(dark) regions show maxima(minima) of 
\J(x,y). Superposed on top are arrows showing the two- 
dimensional vector field J(x,y) where the arrow lengths are 
scaled to \J\. The rotations are produced by inversion of drift 
currents produced by opposing ratchet potentials in a given 
direction. The current density J scales with the asymmetry 
parameter "a" and the square of the correlation time r c so for 
white noise or symmetric potentials, there are no rotations. 



Table I shows the three general cases of multi- 
dimensional flows which may be created by appropriate 
choice of the potential. Figures 4-7 show contour plots of 
different potentials and their corresponding current den- 
sities J. The first case is shown in Fig0| which shows 
that for non-separable potentials one gets in general ro- 
tational vortices. This demonstrates the breaking of ro- 
tational symmetry of the system by construction of an 
appropriate ratchet potential. Notice that the total cir- 
culation along the boundary is zero, which means that 
although there are local circulation patterns, the global 
average is zero. This is a consequence of periodic bound- 
ary conditions that we use. Global rotations will need a 
net rotation along the boundary, which is a different set 
of boundary conditions (similar to that in Fig. than 
we are analyzing here. The absence of any global circu- 
lation causes the current circulation patterns to come in 
vortex-antivortex pairs. 




0.2 0.4 0.6 0.8 1 

x 

FIG. 5: Same as Fig. 21 with the potential replaced with a 
separable one V(x,y) = sin(x) — sin(2x) + sin(y) — sin(2y), 
which leads to laminar flow, caused by two independent de- 
coupled ratchets in the X and Y directions. 




X 



FIG. 6: Combination of drift and laminar flow in coupled 
potential U(x, y) = [sin(x)sin(y) — sin(2x)sin(2y)\ +0.2(:r + 
y). The term in x + y is periodically repeated outside the 
interval [0,1]. 

FigEl shows how in a separable potential one gets lam- 
inar flow. By adjusting the relative magnitude of the 
potential terms involving the x and y coordinates one 
can change the angle the laminar flow makes with the X 
axis. Thus we can obtain the breaking of reflection sym- 
metry about the two axes x and y, with this flow pattern. 
Notice that although there are local swings in the current 
density, there are no vortex-antivortex pairs, since only 
the relative magnitudes (but not the signs) of J x (y) and 
J y {x) change. 

For a combination of separable and non-separable po- 



1 

0.8 
0.6 

z 

0.4 
0.2 




/ 



1 



0.5, 



X 



0.5 

y 




FIG. 7: 3-D flow patterns in the potential \J(x,y,z) — A [sin x sin y sin z — a sin (2x) sin (2y) sin (2z)] for A = 1, a = 1. The 
3-D potential leads to a rich structure of loops and vortices (left). In the panels to the right, components J x and J y of the 
current are plotted at various z slices, corresponding to heights z/2n = (a) 0.2 (top left), (b) 0.4 (top right), (c) 0.6 (bottom 
left) and (d) 0.8 (bottom right). Note that the sizes of the arrows show the relative magnitudes of the J x ,J y currents within 
each panel, and should not be compared for varying heights (such a comparison is shown in the 3-D plot on the left). 



tentials as shown in Fig. |^1 one gets equivalent combi- 
nations of patterns i.e. one can generate net drifts along 
with local rotations. Fig. [7| shows a flow pattern for a 3 
dimensional potential, shown as a 3-D plot, and as 2-D 
plots at different z slices. Infinitely many combinations of 
such potentials can be generated to break different kinds 
of symmetries in higher dimensions as shown in Figs. (4- 
7). As an example, in 3-D one could have a 1-D ratchet 
along the z-dircction completely decoupled from a rotat- 
ing 2-D ratchet in the x-y plane. This would give rise to 
a helix with its axis along the z-direction. 

Finally, we note that we have used Dirichlet boundary 
conditions with specified P(x, 0) and P(0, y) functions 
(assumed to be constant), in order to get our flow pat- 
terns. We could use Neumann boundary conditions as 
well, by rewriting P(0,y) in J x (y) in terms of the nor- 
mal derivative dP(x,y)/dx\ x =o- This is done using the 
definition of J from Eq. [5] and in Eq. |2HI setting = x: 



P(o,y) 



1 e,(o,„, ap <-«" 



x(y) 



Ox 



(33) 



x=0 



where 



x(y) = yVx{x,y) 
- e x (o,y) 



dO x {x,y) 
dx 

[1 _ e -0x(Ar,y)] 



x=0 



(34) 



and an analogous equation for P(x, 0) in terms of 
dP(x,y)/dy\ 



VII. CONCLUSIONS 

We have shown that a Fokker-Planck Equation may 
be derived for weak colored noise for a bounded poten- 
tial in multiple dimensions. Using this we have demon- 
strated that an asymmetric periodic potential with a col- 
ored noise in n dimensions will necessarily lead to break- 
ing of many types of symmetries of the particle motion. 
Specific examples have been solved for laminar flow and 
laminar flow combined with rotations in two and three 
dimensions. Such symmetry-breaking in higher dimen- 
sions should also be readily generalizable to other types 
of noise statistics, such as a non-Gaussian probability 
distribution, and other types of ratchet potentials. This 
could include time-dependent ( "flashing" ) ratchet poten- 
tials [2l| , as well as discrete versions such as multidimen- 
sional Parrondo's games [2(| . 
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IX. APPENDIX (A) : MATHEMATICAL 
DETAILS 

(i) Lemma: 5N/5f k (t) = 
Proof: 



SN 
SMt) 



5 h{t) 



(35) 



From the normalization condition J P[f\T>f = 1 and Eq. 
© we get 



N 



= / Vfexp 



Hence we get 
SN 



\J J dsds'Y^Kij{s- s')fi(s)fj{s') 

(36) ■ 



6fk{t) 



~N 2 J £>/exp(...) 



dsds' 



x ^S ik S{t - s)/^') + S jk S(t - s')fi(s) 

= W Z>/exp (...)£ J ds(Ms)/2) 

x [K ki (t-s)+K ik (s-t)] 

= (N/2) J2j ds [KmQ - s ) + R ik( s - *)] 



x y vfp[f}fi(s) 



(37) 



where the second integral on the last line can be rewritten 
as (fi), which is zero since fi is a noise. This completes 
the proof. 

As a corollary to the above, we obtain the following 
two equations: 



SP[f] 
SfkQk) 



= - E / dsK ^ - s )fi( s ) p lf\ ( 38 ) 



+ E / / dsds'K tk {t ~ s)!^^ 1 - s')f 3 {s')P[f\ 

(39) 

(ii) Lemma: Jds'K u (t' ~ s')Ci([s - s']/rf) = TfS a S(t - 
s')/A =fds'K a (s - s')Ci([t' ~ s']/t9). 

This equation establishes that if is a diagonal matrix, 
whose inverse gives C . Differentiating the normalization 
equation 1 = J VfP[f], we get: 



Sf k (t)5Mt> 
S 2 P[f] 



VfP[f] 



Vf 



-K kl (t-t') 



dsds' 



Kikit-sWtf-s'WitfW)) 

(40) 



Rewriting K kl {t - t') as YhS dsK *k(t - s)S a S(s - t'), 
and using (/<(«)/,•(«')> = ^ (A/rf) Ci([« - s']/rf) leads 
immediately to the above proof. 



X. APPENDIX (B) : CONVERGENCE ISSUES 

It is important to establish the validity of various ex- 
pansions that we do in the correlation time at various 
stages of the derivation of Eq. ^| In particular, the 
current distribution may not necessarily be an analytic 
function of r c for an arbitrary noise-source, in which case 
any perturbative expansion in r c yields results that are 
wrong. The non-analyticity can arise from three possible 
sources: (i) the potential itself may be non-analytic, (ii) 
the noise statistics has a finite support, (iii) the pertur- 
bative expansion may be of a non-analytic function of 
r c . 

Case(i) A piecewise linear sawtooth potential has in- 
finite derivatives at the kinks, leading to potential di- 
vergent terms in the current. This can be avoided in 
principle by restricting our arguments to a smoothened 
function. 

Case(ii) The derivation of the Fokker-Planck structure 
itself depends crucially on the statistics of the noise. If 
the noise distribution has a finite support so that arbi- 
trarily large noise amplitudes are excluded from consider- 
ation, then any computation of the current density along 
the lines we prescribed would be totally wrong. For ex- 
ample, if the height of the potential barrier is larger than 
the maximum allowed noise amplitude, then there will 
be no current, contrary to what an injudicious applica- 
tion of the formalism will yield. In our analyses, we have 
assumed a Gaussian distribution function for the noise 
(Eq. [21 • This has an infinite support and thereby avoids 
such non-analyticities j24[. However, for a discrete noise 
process such as dichotomous noise, care must be exer- 
cised in obtaining the Fokker-Planck description. Often, 
an additional white noise source is included with the ex- 
plicit purpose of handling such non-analyticities. We do 
not need such sources since our probability distribution 
is Gaussian. 

Case(iii) The functional calculus as outlined by Fox 
and extended by us to several dimensions is non- 
perturbative in the correlation time. As argued by Fox 
in Ref. the prescription leads to currents that are 
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uniformly convergent for r c = 0. The non-perturbative 
description leads us to Eq. that involves the exponen- 
tial of integrals of matrix elements of M. We finally 
performed an explicit evaluation of the matrix elements 
in Eq. ^2f° r small r c . This does involve a perturbative 
expansion in r c , but of an exponential function, analytic 



in t c . The nth term of the expansion is proportional to 
Hu{t c /t 7 )™ jn\. For well-behaved correlation functions 
such as an exponential or a Gaussian, this term tends to 
zero rapidly as n increases to infinity, provided t c < r 7 , 
as we have assumed. 
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